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METHOD AND PROGRAM FOR DETERMINING INTERSECTION POINT 



OF TRIANGLE WITH LINE SEGMENT 



BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a method and a program for 
determining an intersection point of a triangle with a line segment. 

Description r of the Related Art 

Substantial data that includes shape data and physical 
property data integrated with each other can be stored with a small 
memory capacity. Thereby, data regarding a shape, a structure, a 
physical property and a history of an object can be managed in an 
integrated fashion, and data regarding a series of processes such 
as designing, machining, assembling, testing, and evaluating can be 
managed by the same data. Patent Literature 1 discloses a substantial 
data storing method by which a CAD and a simulation can be integrated. 
[Patent Literature 1] Japanese Laid-Open Patent Publication No. 
2002-230054 . 

According to Patent Literature 1, as shown in FIG. 1, the method 
of storing substantial data configured from shape data and physical 
property data integrated with each other includes an external data 
input step (A), oct-tree division step (B) , and a cell data storing 
step (C) . In the external data input step (A) , the external data 12 
including the boundary data of the object is input to a computer or 
the like having a stored program for the method of the Patent 
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Literature 1. This boundary data is obtained in an external data 
obtaining step SI. In the oct-tree division step (B) , oct-tree 
division divides the external data 12 into rectangular parallelepiped 
cell of which boundary surfaces are perpendicular to each other. In 
5 the cell data storing step (C) , various physical property values are 
stored for each cell. 

According to the storing method of Patent Literature 1, the 
external data including the shape data of the object is divided into 
rectangular parallelepiped cells of which boundary surfaces are 
10 perpendicular to each other, and the various physical property values 
are stored for each cell. Each divided cell is classified into an 
internal cell or a boundary cell that includes the boundary surface. 
The internal cell contains one kind of physical property values as 
attribute, and the boundary cell contains two kinds of physical 
15 property values corresponding to the inside and the outside of the 
object, respectively. 

Hereinbelow, the data used in the above method is referred to 
as "V-CAD data", and designing and simulation using this data is 
referred to as "volume-CAD" or "V-CAD". In FIG. 1, the reference 
20 numeral 14 designates V-CAD data. 

FIG. 2 shows a flow of creating V-CAD data from external data 
according to the method of Patent Literature 1. As shown in FIG. 2, 
in STEP-1, a voxel space is defined for S-CAD data or a triangle patch 
data, and in STEP-2, intersection points of a geometric shape with 
25 the ridges of each cell are determined by calculation (that is, cell 
cut points are determined) . In STEP-3, an outer circumferential loop 
of the cell internal surface is formed on the basis of the cut point 



data for each cell . In STEP-4 , the inside of the outer circumferential 
loop are divided into triangles by referring to the original geometric 
shape. In STEP-5, the medium values of the cell are set. Thereby, 
the V-CAD data is formed. 

In the method of Patent Literature 1, the intersection points 
of the triangle with the ridges of each cell need be . repeatedly 
calculated a number of times. This triangle defines a geometric shape 
expressed by the external data. 

Non-Patent Literatures 1 through 3 disclose conventionally 
known methods of calculating the intersection point of the triangle 
with the line segment. 
[Non-Patent Literature 1] 

Moller's method (Journal of graphics tools, 2(1): P21-28 

1997) . 

[Non-Patent Literature 2] 

Badouel's method (Graphic Gems, Academic Press P21-28 1990. 
[Non-Patent Literature 3] 

Computational Geometry in C: Joseph 0 1 Rourke/Cambridge 
University Press P239) . 

Various methods are used for calculating the intersection 
point of a triangle with a line segment. Many of these methods 
calculate the intersection of an infinite plane with an infinite line 
to determine whether or not the intersection point i,s positioned 
inside the triangle of the closed region on the plane . The calculation 
by these methods results in a high cost, and in the future in the 
field of the V-CAD system, it will be necessary to improve the 
calculation method in which calculation is repeatedly performed a 
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number of times. 

According to the Moller's method of Non-Patent Literature 1, 
a calculation cost can become relatively low by use of the vector 
calculation expressed by the following Equation (1) of Formula 2. 
[Formula 2] 



Lin 



(1) 



Pi/ P2/ P3/ PLin: the triangle vertexes, the starting point of 
the line segment 
10 V Dir : a direction vector of the line segment 

t: a parameter in the line segment 
u and v: parameters in the triangle 

However, in this method, even when calculating intersection 
points of plural different line segments with the same triangle, it 
15 is necessary to work out the equation repeatedly for each line segment . 
Accordingly, this method produces no merit for the intersection 
calculation on the same triangle, and cannot improve a calculation 
speed. 

According to the Badouel's method of Non-Patent Literature 2, 
20 an intersection point is calculated by using a standard point P Pln 
of a plane, a normal vector V Nrm , and a standard point P Lin of a line. 

The distance between the standard point of the line and the 
plane is calculated by dst - V Nrm • V Pnt (V Nrm is a unit vector) , so that 
the coordinates of the intersection point are calculated by P Int = 
25 P Lin + V Lin • dst/ ( 

v Lin rra / • 
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Subsequently, the respective triangle vertexes and the 
intersection point of the plane with the line are projected on any 
of coordinate planes (xy plane, yz plane and zx plane), and a and 
B in the equation (3) of Formula 3 are calculated on this projection 
5 plane. 

In the case of the xy projection plane in FIG. 4, a and j8 are 
calculated by the inverse transformation of the linear transformation 
of the equation (4) of Formula 3 with Pi(xi, y lt z x ) , P 2 (x 2 , y 2 , z 2 ) 
and P 3 (x 3 , y 3 , z 3 ) and P Int (x, y, z) , and with the z coordinates of each 
10 vertex being omitted. 
[Formula 3] 
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Vpnt^PpinPLin ' * * (2) 

P^P£t= 0(F&'2 + FP r fc , 3 ' ' * (3) 

wV» *-4p) ■ * ' (4) 

There is following another inside/outside determination 
15 method of determining whether or not the intersection point is inside 
the triangle while this inside/outside determination method 
determines the intersection point of an infinite plane with a line 
by the same calculation as that of the Badouel ? s method. 

According to an angle calculation method of Non-Patent 
Literature 3, in FIGS . 5A and 5B, inner product is calculated for 
each two of vectors extending from the intersection point to the 
respective vertexes. Then, it is determined whether or not the 
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intersection point is inside the triangle on the basis of the 
calculated angle sum of 0 or 360 degrees. That is, in the case of 
FIG. 5A, the equation (5) can be established, and in the case of FIG. 
5B, the equation (6) can be established. 
5 0i 2 +0 23 + 0 31 =2PI ... (5) 

012+023+031 = 0 ... (6) 

Alternatively, as shown in FIGS. 6A and 6B, outer product is 
calculated for each two of the vectors from the intersection point 
to the respective vertexes by the equation (7) of Formula 4. 
10 [Formula 4] 

VNrml2 — Vl X Vfc 

VNrm23 = V 2 X V 3 • • • (j ) 

VNrm31 = V 3 X Vl 

In the case of FIG. 6A, V Nrm i 2 , V Nrm2 3 and V Nrm31 extend in the same 
direction, and in the case of FIG. 6B, VNrm31 extends in the opposite 
15 direction to V Nrml2 and V Nrm23 . 

However, calculation of the outer product will apparently 
increase calculation load. 

As described above, according to the conventional method of 
calculating the intersection point of a triangle with a line segment, 
20 when calculating the intersection point of the same triangle with 
plural different line segments, the calculation need be performed 
again and again for each intersection point. As a result, entire 
calculation time becomes long, and it becomes difficult to realize 
high speed calculation. 

25 



SUMMARY OF THE INVENTION 



The present invention was made to solve the above problems. 
That is, it is an object of the present invention to provide 
intersection point calculation method and program that can reduce 
the calculation time and achieve high speed calculation at the time 
of calculating intersection points of the same triangle with plural 
different line segments. 

According to the present invention, there is provided a method 
for calculating an intersection point of a triangle with a line segment, 
the triangle having three vertexes P 0 , Pi, P 2 , and the line segment 
having both ends A and B, comprising: 

a coordinate system setting step (A) of setting a coordinate- 
system R in which P 0 is a coordinate origin, P 0 Pi conforms to a first 
axis (U axis) to have a unit length, P0P2 conforms to a second axis 
(V axis) to have a unit length, and P 0 Pi x P 0 P 2 is a unit vector conforming 
to a third axis (N axis) , 

a transforming matrix calculation step (B) of calculating a 
transforming matrix M that transforms coordinate values in an ordinary 
coordinate system into coordinate values in the coordinate system 
R; 

a line segment end coordinate calculation step (C) of calculating 
u-, v- and n-coordinate values of the both ends A and B of the line 
segment; 

an intersection determination step (D) of determining whether or not 
the line segment intersects with the triangle, on the basis of the 
u ~/ v- and n-coordinate values of the both ends A and B; 
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an intersection point coordinate calculation step (E) of 
calculating u-, and v-coordinate values of the intersection point 
when the line segment intersects with the triangle; and 

an intersection point determination step (F) of determining 
5 whether or not the intersection point is positioned inside the 
triangle, on the basis of the u-, and v-coordinate values of the 
intersection point. 

Further, according to the present invention, there is provided 
a program causing a computer to calculate an intersection point of 
10 a triangle with a line segment, the triangle having three vertexes 
Po, Pi/ P2/ and the line segment having both ends A and B, 

the program causing a computer to perform: 

a coordinate system setting step (A) of setting a coordinate 
system R in which P 0 is a coordinate origin, P 0 Pi conforms to a first 
15 axis (U axis) to have a unit length, P0P2 conforms to a second axis 
(V axis) to have a unit length, and P0P1 x P0P2 is a unit vector 
conforming to a third axis (N axis), 

a transforming matrix calculation step (B) of calculating a 
transforming matrix M that transforms coordinate values in an ordinary 
20 coordinate system into coordinate values in the coordinate system 
R; 

a line segment end coordinate calculation step (C) of calculating 
u-, v- and n-coordinate values of the both ends A and B of the line 
segment; 

25 an intersection determination step (D) of determining whether or not 
the line segment intersects with the triangle, on the basis of the 
u-, v- and n-coordinate values of the both ends A and B; 
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an intersection point coordinate calculation step (E) of 
calculating u-, and v-coordinate values of the intersection point 
when the line segment intersects with the triangle; and 

an intersection point determination step (F) of determining 
5 whether or not the intersection point is positioned inside the 
triangle, on the basis of the u-, and v-coordinate values of the 
intersection point. 

With the above method or program, it is possible to calculate 
the u-, v- and n-coordinate values of the both ends A and B of the 
10 line segment by only performing the line segment end coordinate 
calculation step (C) on the plural different line segments because 
the transforming matrix corresponding to the same triangle is always 
same when intersection points of the same triangle with the plural 
different line segments are calculated. Accordingly, the present 
15 invention can reduce the number of calculations compared to the 
conventional method so that the calculation time can be reduced, and 
high speed calculation can be achieved. 

According to a preferred embodiment of the present invention, 
the transforming matrix M calculated in the transforming matrix 
20 calculation step (B) is an affine transforming matrix. 

The affine transforming matrix corresponding to the triangle 
can easily be calculated in a short time. 

Furthermore, preferably, in the line segment end coordinate 
calculation step (C) , the u-, v- and n-coordinate values are 
25 calculated by the equation (8) : 

P 0VN = MP xyz . . . (8) . 

By using the equation (8), it is possible to easily calculate 
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the u-, v- and n-coordinate values of the both ends A and B of the 
plural different line segments in a short time. 

Preferably, the intersection determination step (D) 
including: 

5 (1) determining that the intersection point does not exist when 

one of the both ends has the absolute n-coordinate value not smaller 
than a length of the line segment, and/or 

(2) determining that the intersection point exists when the 
n-coordinate values of the both ends have different signs. 
10 Thereby, it is possible to easily determine whether or not the 
intersection point exists, in a short time. 

Preferably, the intersection determining step (D) includes 
determining that the intersection point is positioned outside the 
triangle (1) when each of the both ends has the u-coordinate value 
15 not larger than "0 M or not smaller than "1", (2) when each of the 
both ends has v-coordinate value not larger than "0", or (3) when 
each of the both ends has the u- and v-coordinate values of which 
sum is not smaller than "1". 

Thereby, it is possible to easily determine whether the 
20 intersection point is inside or outside the triangle, in a short time. 
Further, it is possible to exclude the line segment that does not 
apparently have the intersection point inside the triangle, omitting 
superfluous calculations . 

Preferably, the intersection point coordinate calculating 
25 step (E) includes internally dividing the line segment by the n- 
coordinate values of the both ends A and B, determining the internally 
divided point as the intersection point, and calculating u- and 
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v-coordinate values of the intersection point. 

Thereby, it is possible to easily calculate the u- and v- 
coordinate values of the intersection point in a short time. 

Furthermore, preferably, the intersection point determination 
5 step (F) includes determining that the intersection point is inside 
the triangle when the u- and v-coordinate values of the intersection 
point satisfy three inequalities (13) of Formula 1: 
[Formula 1] 

• 0^ u (range in a Ved direction) 

• v ^ 1 (range in a Vec2 direction) • • ■ (13) 

• 0 ^ u+v ^ l (inside a line Pi R2) 

10 Thereby, it is possible to easily determine whether or not the 

intersection point is inside or outside the triangle, in a short time. 
Further, before the intersection point is actually calculated, the 
line segment that does not have the intersection point inside the 
triangle can be excluded, omitting superfluous calculations. 

15 Other object and advantageous features of the present 

invention will become apparent in the following description with 
reference to the attached drawings. 

BRIEF DESCRIPTION OF DRAWINGS 

20 

FIG. 1 is a flowchart showing a substantial data storing method 

• * 

of a prior application; 

FIG. 2 is a flowchart showing another substantial data storing 
method of a prior application; 
25 FIG. 3 is a schematic illustration showing a conventional 
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method; 

FIG. 4 is a schematic illustration showing another 
conventional method; 

FIGS. 5A and 5B are schematic illustrations showing another 
5 conventional method; 

FIGS. 6A and 6B are schematic illustrations showing another 
conventional method; 

FIG. 7 shows a configuration of an apparatus that performs the 
method of the present invention; 
10 FIG. 8 is a flowchart showing the method and program of the 

present invention; 

FIG. 9 is a schematic illustration showing a coordinate system 
R according to the present invention; 

FIG. 10 is a flowchart concretely showing processes according 
15 to the present invention; 

FIG. 11 is a flowchart showing a first embodied example 
according to the present invention; 

FIG. 12 is an image of input triangles; 

FIG. 13 is an image showing V-CAD data constituted by a cell 
20 and the internal surface of the cell; 

FIG. 14 is an image showing cell cut points produced by 
calculating an intersection point between input data and cell ridge; 

FIG. 15 is a flowchart showing a second embodied example 
according to the present invention; and 
25 FIG. 16 is an image showing V-CAD data after shape modification 

process . 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following, a preferred embodiment of the present 
5 invention will be described with reference to the drawings. 

FIG. 7 is a configuration of an apparatus for implementing the 
method of the present invention. As shown in FIG. 7, this apparatus 
10 includes an input device 2, an external storage device 3, an 
internal storage device 4, a central processing unit 5 and an output 
10 device 6. 

The input device 2 is a keyboard, for example, and inputs 
external data 12 including shape data of an object 1. The external 
storage device 3 is a hard disk, a floppy disk, a magnetic tape, a 
compact disk or the like, and stores substantial data of shape data 

15 and physical property data integrated with each other, and a storage 
operation program for the substantial data. The internal storage 
device 4 is a RAM or a ROM, for example, and stores the operational 
information. The central processing unit 5 (CPU) mainly processes 
the operation, input and output, and performs the storage operation 

20 program in cooperation with the internal storage device 4 . The output 
device 6 includes a display and a printer, for example, and outputs 
the stored substantial data and the operational result of the storage 
operation program. 

FIG. 8 is a flowchart of the intersection point calculation 

25 method and program of the present invention. As shown in FIG. 8, the 
intersection point calculation method and program include a 
coordinate system setting step (A) , a transforming matrix calculation 
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step (B),a line segment end coordinate calculation step (C) , an 
intersection determination step (D) , an intersection point coordinate 
calculation step (E) , and an intersection point determination step 
(F) . 

5 The intersection point calculation method of the present 

invention is a method of calculating the intersection point of a 
triangle with a line segment, the triangle having vertexes P 0/ Pi and 
P 2 , and the line segment having both ends A and B. The intersection 
point calculation program of the present invention causes a computer 
10 to calculate the intersection point of a triangle with a line segment, 
the triangle having vertexes P 0 , Pi and P 2 , and the line segment having 
both ends A and B. 

In the coordinate system setting step (A) , a coordinate system 
R is set. In the coordinate system R, P 0 is a coordinate origin, P 0 Pi 
15 conforms to a first axis (U axis) to have a unit length, P0P2 conforms 
to a second axis (V axis) to have a unit length, and P 0 Pi x P0P2 is 
a unit vector that conforms to a third axis (N axis). 

In the transforming matrix calculation step (B) , a 
transforming matrix M is calculated. The transforming matrix M 
20 transforms coordinate values in a normal coordinate system into 
coordinate values in the coordinate system R. 

In the line segment end coordinate calculation step (C) , u-, 
v- and n-coordinate values of the both ends A and B are calculated. 
In the intersection determination step (D), it is determined whether 
25 or not the line segment intersects the triangle, on the basis of u-, 
v- and n-coordinate values of the both ends A and B. 

In the intersection point coordinate calculation step (E) , u-, 
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and v-coordinate values of the intersection point are calculated when 
the line segment intersects the triangle. 

In the intersection point determination step (F), it is 
determined whether or not the intersection point is positioned inside 
5 the triangle, on the basis of u-, and v-coordinate values of the 
intersection point. 

Next, the present invention will be described in detail. 
1. In development of the volume CAD, the intersection point 
of the line segment with a triangle patch is repeatedly calculated. 

10 This triangle patch is often used for an input surface shape. In the 
calculation of the intersection point of the triangle with the line 
segment (cell ridge) , an intersection point of a plane with a line 
is calculated, and it is then determined whether or not the 
intersection point is positioned inside an effective region on the 

15 line, and whether or not the intersection point is positioned inside 
.an effective region on the plane (this determination is referred to 
as an inside/outside determination) . Thus, a vector operation is 
performed many times, so that a calculation load becomes large. For 
example, for one kind of input triangle patch group, when the input 

20 triangle patch group is converted into volume data such as V-CAD (e.g., 
data expressing a three-dimensional entire extent, or data for CT, 
MRI or volume rendering) , the division number and the resolution are 
changed, or the input data becomes lattice data having a hierarchical 
structure such as Oct-tree data or multiplex voxel data. In such a 

25 case, intersection points of one surface with edges of the lattice 
are calculated simultaneously or repeatedly one by one, resulting 
in a large amount of calculation. Accordingly, the conventional 
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calculation method is not efficient, and a calculation method that 
substantially improves efficiency is desired. In the future, not 
only volume CAD but also other volume data are widely used because 
of data abundance, and for use of the volume data, a combination of 
5 a triangle patch and volume data is desired because of easy process 
and control. The present invention provides important technique for 
the combination of the triangle patch and the volume data. 

According to the present invention, the input data is 
transformed by an af fine transforming matrix which is generally used 

10 for displaying a shape on a computer. Thereby, the intersection point 
calculation and the inside/outside determination can be performed 
stably at a high speed. Particularly, in the case of repeatedly 
calculating the intersection point for the same triangle, the already 
created triangle-peculiar matrix is re-used to greatly reduce a 

15 calculation cost. For example, a mechanical component is often 
expressed by planes, so that a large triangle is used for input surface 
shape, and an intersection point of the same triangle with a cell 
ridge is calculated a number of times. In such a case, the present 
invention is outstandingly advantageous. 

20 

2. According to the present invention, as shown in FIG . 9, 

the two sides of the triangle are set to conform to first and second 

axes, respectively, and the normal of the triangle is set to conform 

« « 

to a third axis. The triangle vertex shared by the first and second 
25 axis is set to be the origin of the coordinate axes. In this manner, 
the space coordinate system R can be established. 

To convert the orthonormal coordinate system into the 
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above-described space coordinate system, an affine transforming 
matrix which is used for displaying an image on a computer, for example 
is calculated so as to relate the orthonormal coordinate system with 
the triangle. 

The coordinates of the both ends of the line segment for which 
the intersection point is calculated are transformed into the 
coordinates in the space coordinate system R. Then, on the basis of 
the values of the transformed coordinates, it is determined whether 
or not the line segment exists, and whether or not the intersection 
point is positioned inside the triangle. 

The use of the above method produces the following advantages: 

(1) The intersection point calculation and the inside/outside 
determination can be performed by one step without increasing a 
storage capacity, enabling high speed calculation while the 
intersection point calculation and the inside/outside determination 
can be performed by two steps in the conventional method. 

(2) When calculating the intersection points of plural line 
segments with the same triangle, a previously prepared geometric 
property is used, reducing a calculation cost. 

(3) Since a simple algorism for the coordinate transformation 
and the inner product is used, the method can easily be embodied. 

(4) The data obtained by the coordinate transformation 

includes the signed distance so that this signed distance can be used 

< 

for calculating the distance from the original curved surface as the 
attribute data, and for confirming the front and back directions of 
the surface. 
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3. FIG. 10 is a flowchart concretely showing processes 
according to the present invention. As shown in FIG. 10, the 
three-dimensional coordinate data of the triangle having the vertexes 
Po, Pi and P 2 is input. Then, in the above-described coordinate system 
5 setting step (A) , the coordinate system R is produced in which P 0 is 
the coordinate origin, the line segment P0P1 conforms to the first 
axis (U axis) to have a unit length, P0P2 conforms to the second axis 
(V axis) to have a unit length, and P0P1 x P0P2 is a unit vector that 
conforms to the third axis (N axis) 
10 Next, at the step SI, a matrix corresponding to the coordinate 

system R is calculated on the basis of the triangle. This step SI 
corresponds to the above-described transforming matrix calculation 
step (B) . 

Subsequently, at the step S2, the u and v coordinates of the 
15 both ends A and B of the line segment, and the signed distance of 
the both ends are calculated by the affine matrix operation. This 
step S2 corresponds to the above-described line segment end coordinate 
calculation step (C) . 

Next, at the step S3, it is determined whether or not the line 
20 segment intersects with the triangle, on the basis of the signed 
distance and the u and v coordinates. This step S3 corresponds to 
the above-described intersection determination step (D) . 
Next, at the step S4, the intersection point is calculated. This step 
S4 corresponds to the above-described intersection point coordinate 
25 calculation step (E) . 

Then, at the step S5, it is determined whether or not the 
intersection point is positioned inside the triangle, on the basis 
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of the u and v coordinates of the intersection point. This step S5 
corresponds to the above-described intersection point determination 
step (F) . 

5 4.1 The calculating procedure will be described in detailed. 

4.1.1 The coordinate system based on the triangle, and the matrix 
corresponding to the triangle. 

As shown in FIG. 9, on the assumption that the vertexes of the 
triangle are the points P 0 , Pi and P 2 , the coordinate system R is 

10 produced such that in the coordinate system R, P 0 is the coordinate 
origin, P 0 Pi conforms to the first axis (U axis) to have a unit length, 
P0P2 conforms to the second axis (V axis) to have a unit length, and 
P 0 Pi x P0P2 is a unit vector that conforms to the third axis (N axis) . 
Assuming that M is the matrix that transforms the ordinary orthonormal 

15 coordinates into the coordinates in the coordinate system R, the 
following equation (8) : 

Puvn = MPxyz . . . ( 8 ) 
produces UV values on the coordinate surface, and the signed distance 
from the triangle that are calculated as u, v and n coordinates, which 

20 are determined from the triangle. 

Thereby, it is possible to calculate or determine position 
relation between the triangle and an arbitrary point. For example, 
it is possible to calculate the signed distance between the triangle 
and the arbitrary point, or determine whether or not the point can 

25 be orthogonally projected on the triangle. Furthermore, when the 
intersection point is calculated a number of times for the same 
triangle, the matrix can be used without modification, so that 
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calculation load can be reduced. As a result, the intersection point 
can be calculated at the higher speed. 

Since this matrix includes functions of scaling, rotation, 
coordinate axis change, and parallel movement, transforming matrix 
5 of 4 x 4 is used. This transforming matrix is usually used for 
displaying computer graphics, and is generally referred to as an 
af f ine transformation . 

If respective axis direction scalings are expressed by xScale, 
yScale and zScale, and the parallel movement is expressed by (xMove, 
10 yMove, zMove) , the scaling matrix and the parallel movement matrix 
can be respectively expressed by the equations (9) and (10) of Formula 
5. 

The unit vectors in the respective axis directions in the uvn 
coordinate system are Vecl, Vec2 and VecN, respectively. If the 
15 components of Vecl are (vxl, vyl, vzl) , the components of Vec2 (vx2, 
vy2, vz2), and the components of Vec3 (vxn, vyn, vzn) , the rotation 
matrix and the coordinate change matrix can be expressed by the 
equation (11) of Formula 5. 

The transforming matrix is formed as the equation (12) of 
20 Formula 5 by unifying the scaling matrix, the parallel movement matrix 
and the rotation and coordinate axis change matrix. 
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[Formula 5] 
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4.1.2 The coordinate transformation of the both ends of the line 
segment (calculation of signed distance and uv coordinates) . 

The affine transformation is performed on the start point and 
end point of the line segment by using the above-mentioned 
transforming matrix. In the affine transforming matrix of 4 x 4, each 
fourth column of the first to third lines is a movement amount, so 
that the multiplication target is always "1" . In the matrix, the first 
to third columns of the fourth line is "0", the fourth column of the 
forth line is "1", and the multiplication target of the fourth column 
is always "1", so that the calculation of the fourth line always 
becomes "1". Accordingly, if instead of usual matrix calculation, 
the affine transformation calculation is performed, the calculation 
is performed for only the matrix components ranging from the firth 
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line to third line and from first column to the third column, so that 
the necessary number of the calculation becomes only nine . As a result , 
the high speed calculation can be realized. 

4.1.3 Determination of the existence of the intersection point on 
5 the basis of coordinates of two transformed points. 

If one of the two transformed points has an absolute n-coordinate 
value larger than the length of the line segment, the intersection 
point does not exist. 

If the n coordinates of the two transformed points have the 
10 same sign, the intersection point does not exist, and if the n 
coordinates of the two transformed points have different signs, the 
intersection point exists. 

4.1.4 Determination of whether the possible intersection point is 
inside the triangle, on the basis of u and v coordinates of two 

15 transformed points 

If both of the u-coordinate values of the two transformed points 
are equal to or smaller than "0", or equal to or larger than f, l", 
the intersection point is outside the triangle. 

If both of the v-coordinate values of the two transformed points are 
20 equal to or smaller than "0.0", the intersection point is outside 
the triangle. 

If the sum of the u-coordinate value and v-coordinate value 
of each of the two transformed points becomes equal to or larger than 
"1.0", the intersection point is outside. the triangle. 
25 These requirements can exclude the line segment that does not 

apparently have the intersection point before actually calculating 
coordinate values of the intersection point. 
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4.1.5 u- and v-coordinates are calculated on the assumption that the 
intersection point is the point internally divided by the n-coordinate 
values of the both ends of the line segment. 

4.1.6 It is determined whether or not the intersection point is 
5 positioned inside the triangle, by using the following requirement 

for the u and v coordinates. 

The intersection point is inside the triangle if the following 
requirement of the bottom one of inequalities (13) of Formula 6 is 
satisfied. 
10 [Formula 6] 

• 0 ^ u ^ l (range in a Ved direction) 

• 0 ^ v ^ l (range in a Vec2 direction) • • • (13) 

• 0 ^ u+v ^ l (inside a line P1R2) 

4.2 Calculation cost and calculation result 

The calculation of the present invention is compared to that 
15 of the Moller's method. 

In the Moller's method, each vector product in the equation 
(14) of Formula 7 is used two times in the equation (15) of Formula 
7. Accordingly, these vector products are previously calculated to 
be vectors P and Q, respectively. The division calculation in the 
20 equation (15) is previously performed to obtain the reciprocal. 

In this case, the number of the multiplication calculations 
for obtaining the vector P becomes 6, and the number of the 
multiplication calculations for obtaining the vector Q becomes 6. 
Accordingly, the calculations for obtaining the value of "t" includes 
25 the multiplication for obtaining the vectors P and Q, and the inner 



24 



product of the vector P and another vector, and the inner product 
of the vector Q and another vector, reciprocal calculation, and 
multiplication with the reciprocal. That is, to obtain the value of 
"t", the multiplication is performed 19 times, and division 
calculation is performed once. 

To obtain a value "u", the additional 4 multiplications are 
required, and to obtain a value "v", the additional 4 multiplications 
are required. 

On the other hand, according to the present invention, to obtain 
a signed distance and values of the both end points "u" and "v", 18 
multiplications are required. If the signed distance has the same 
sign at the both end points, the value "t" is smaller than "0" or 
larger than "1", and the intersection point does not exists in the 
effective range of the line segment. Accordingly, for the value "t", 
the present invention can save one multiplication and one division 
calculation . 

According to the present invention, to obtain the values "u" and "v", 
the calculation of the equation (16) of Formula 7 is required, so 
that to obtain the value "u M , one division calculation and two 
multiplications are required, and to obtain the value "v", additional 
two multiplications are required. It should be noted that without 
obtaining the values "U" and "V" of each of the start point and the 
end point, it is possible to determine that the line segment does 
not intersects with the triangle if the each of the values "U" and 
"V" of The start point or the end point is smaller than "0", or the 
sum of the values "U" and "V" is equal to or larger than "1" . 
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[Formula 7] 



(PtPu.xPA) ±<y»,*PA) ' ' • (14) 



t ' 




u 




V 





(15) 



t = fabs(N Pl /(N P2 -N Fl )) 

U Int =U FO xt + U Pi x(l-0 • • • (16) 



Accordingly, the intersection calculation is performed in the 
following procedure: 

at Step 1, the existence of the intersection point is determined 
in the effective range of the line segment (t VALUE) ; 

at Step 2, by the value "u", it is determined whether or not 
the intersection point is effective; and 

at Step 3, by the value "v", it is determined whether or not 
the intersection point is effective. 

In the case of this procedure, to determine the existence of 
the intersection point inside the triangle, calculation cost shown 
in Table 1 is required. Thus, the present invention can achieve the 
higher calculation efficiency. 
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[Table 1] 



STEP 


Moller's Method 


This Method 


Stepl tValue 


19Multi+lDivide 


18Multi 


Step2 uValue 


23Multi+lDivide 


20Multi+lDivide 
(l8Multi) 


Step3 vValue 


27Multi+lDivide 


22Multi+lDivide 
(l8Multi) 



4 . 3 Embodied Example 
5 4.3.1 Embodied Example 1 

FIG. 11 shows a procedure from the obtaining of the external 
data to the creating of the V-CAD data. In FIG. 11, before the 
above-described calculation of the intersection point, a process of 
mapping geometrical shape (S-CAD data or triangle patch) to each cell 

10 is performed. By this process, it is possible to reduce the target 
line segments (cell ridge) of the intersection, so that the time 
required for calculation can be reduced, and high speed calculation 
can be realized. If the mapping process is not performed, it becomes 
necessary to search for all line segments, so that the calculation 

15 time increases. For this reason, as in this embodied example, it is 

preferable to perform the mapping process before calculating the 

intersection point. 

FIGS. 12 through 14 are displayed images when the processes 

of FIG. 11 are performed. Specifically, FIG. 12 is the displayed image 

20 of the triangles that are input at STEP T2, and FIG. 13 is the displayed 
« « « 

image of the V-CAD data created at STEP T13. FIG. 14 is an enlarged 
image of the displayed image of FIG. 13, and shows the cut points 
determined by calculating the intersection point between the input 
data and the cell ridge. 
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The calculation method of the present invention enabled the 
calculation speed to be raised by about 30 percent when the benchmark 
was performed. 
4.3.2 Embodied Example 2 
5 FIG. 15 shows a procedure when the method of the present 

invention is applied to an interference check when the shape 
modification process is performed. At STEP U4, the self-interfering 
check is performed by the intersection point calculation (the 
processes of STEP T5 through Til of FIG. 11) before a series of the 
10 processes of FIG. 11 are performed. Thereafter, a series of processes 
of FIG. 11 are performed at STEP U5 so that the modified shape can 
be displayed. FIG. 16 is the displayed image of this modified shape. 
4.3.3 

As described above, when the intersection points of plural 
15 different line segments with the same triangle are calculated, the 
line segment calculation method and program exhibits an excellent 
advantage in that the calculation time can be reduced to achieve high 
speed calculation . 

In the above, the several preferred embodiments of the present 
20 invention are described. However, the scope of the present invention 
is not limited to these embodiments, and the scope of the attached 
claims includes all improvements and equivalents of the present 
invention described in claim. 



